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ABSTRACT 

Shock breakout is the brightest radiative phenomenon in a Type II supernova (SN). Although 
it was predicted to be bright, the direct observation is difficult due to the short duration and X- 
ray/ultraviolet-peaked spectra. First entire observations of the shock breakouts of Type II Plateau SNe 
(SNe IIP) were reported in 2008 by ultraviolet and optical observations by the GALEX satellite and 
supernova legacy survey (SNLS), named SNLS-04D2dc and SNLS-06Dljd. We present multicolor light 
curves of a SN IIP, including the shock breakout and plateau, calculated with a multigroup radiation 
hydrodynamical code STELLA and an evolutionary progenitor model. The synthetic multicolor light 
curves reproduce well the observations of SNLS-04D2dc. This is the first study to reproduce the 
ultraviolet light curve of the shock breakout and the optical light curve of the plateau consistently. 
We conclude that SNLS-04D2dc is the explosion with a canonical explosion energy 1.2 x 10 51 ergs 
and that its progenitor is a star with a zero-age main-sequence mass 20M Q and a presupernova radius 
8OOi?0. The model demonstrates that the peak apparent B-band magnitude of the shock breakout 
would be me ~ 26.4 mag if a SN being identical to SNLS-04D2dc occurs at a redshift z = 1, which can 
be reached by 8m-class telescopes. The result evidences that the shock breakout has a great potential 
to detect SNe IIP at z > 1. 

Subject headings: shock waves — radiative transfer — supernovae: general — supernovae: individual 
(SNLS-04D2dc) — stars: evolution 



1. INTRODUCTION 

In a supernova (SN) explosion, an outward Shockwave 
formed in an inner layer propagates through the stel- 
lar envelope. When the shock emerges from the stellar 
surface, a hot fire ball suddenly appears and flashes in 
soft X-ray or ultraviolet (UV). The flash has a quasi- 
blackbody spectrum (T ~ 10 6 — 10 5 K) and lasts a 
few seconds to ~ 1 days, depending on an ejecta mass 
M e j , e xplosion energy E, and pre supernova radius -R prc sN 
(e.g., iMatzner &: McKeel [1999). The phenomenon is 
called "shoc k breakout" having bee n theoretically pre- 
dicted fe.g. JKlein &: Chevalierlfl978[ ). 

Owing to the short duration and X-ray/UV-peaked 
spectra, it is difficult to observe directly the shock break- 
out. The first detection of shock break out was reported 
for nearby Type II-pecul iar SN 1987A (jCatchpole et al.1 
Il987t fHamuv et aTlll988h . but it was only a detection of 
a rapid decline presumed to be a shock breakout tail. Al- 
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though the detections of shock breakou t tails were occa- 
sionally reported for Type lib SN 1993 J (Richmon d et al.l 
119941 ) and Type lb SN 1999ex (jStritzinger et al.l 12002) . 
there were no observations of shock breakouts from the 
rising part or with the X-ray or UV light. 

First entire observations of the shock breakouts were 
reported in 2008. A first example is an X-ray observation 
of Type lb SN 2008D which fortunately appeared in the 
same galaxy a s SN 2007uy during the observation by the 
Swift satellite (iSoderberg et al.l20 08: Maz zali et al.ll2008l : 
iModiaz et"aTll2009l ; lMalesani et al.ll2009fh Second exam- 
ples are Ty pe II plateau SNe (SNe IIP) SNLS-04D2dc 
(z = 0.185. ISchawinski et al.l [20081: iGezari et al.|[200l 
and SNLS-06Dljd [z = 0.324. ICezari et al.ll2008h . They 
were caught coincidentally i n the UV Deep Imaging Sur- 
veyby the GALEX satellite (jMorrissev et al.l l2005. 2007) 
at the location w here supernova legacy survey (SNLS, 
lAstier et afll2006l ) found SN candidates. 

Since SNe IIP are the most common among core- 
collapse types of SNe (e.g., iMannucci et al.1 120081 ) and 
their shock breakouts are suggeste d to be so bright a s 
to be detected even at z > 1 (e.g., iChugai et al.ll2000T ). 
it is important to develop a way to derive SN properties 
from the shock breakouts. Therefore, we calculate mul- 
ticolor light curves (LCs), including a shock breakout 
and plateau, of SN IIP based on an evolutionary progen- 
itor model with a multigroup radiation hydrodyn amics 
code STELLA (IBlinnikov et al.lH998L [2000. 12006). The 
synthetic LCs are compared with the multicolor obser- 
vations of SN IIP SNLS-04D2dc (R.A. = 10 h 00 m 16.7 s , 
deck = +02°12'18.52" [J2000.0]). 8 We first present a 

8 We focus on SNLS-04D2dc because SNLS-06Dljd has sparse 
UV observations with relatively low signal-to-noise ratio. 
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multicolor LC model reproducing well the shock breakout 
and plateau consistently and constrain SN and progeni- 
tor properties. Furthermore, based on the multicolor LC 
model, we present an apparent B-band light curve of a 
shock breakout of a SN being identical to SNLS-04D2dc 
at z = 1. 

In § [21 the applied models and the radiation hydrody- 
namics calculations are described. In § [3j the multicolor 
LCs of SNLS-04D2dc are compared with the synthetic 
LCs. In § 21 conclusion and discussion are presented. 

2. METHODS & MODELS 

We apply the multigrou p spherical radiation hydro- 
dyna mics code STELLA (jBlinnikov et alj 119981 120001 
2006). STELLA adopts variable Eddington factors, a 
gray transfer of 7-ray from radioactive nuclei, LTE ion- 
ization states, and a multigroup expansion opacity and 
solves the time-dependent equations implicitly for the an- 
gular moments of intens ity averaged over fixed frequency 
bands (for details, see Blinniko v et al.l 120061 and refer- 
ences therein). Multigroup radiative transfer is coupled 
with hydrodynamics, which enables to acquire the spec- 
tral energy distributions (SEDs) self-consistently. The 
color temperature of a SN is estimated from a black-body 
fitting of the SED. In this Letter we adopt 100 frequency 
bins dividing logarithmically from A = 1 A to 5 x 10 4 A; 
the large number of frequency bins allows to describe 
accurately a non-equilibrium continuum radiation. 

A progenitor model is a non-rotating solar-metallicity 
star constructed by a stellar evolution calculation 
(|Umeda fe Nomoto I 120051) . The cal culation includes 
a metallicity-dependent mass loss (jKudritzkil [2000) 
and thus the presupernova model has a self-consistent 
-RpreSN, luminosity, temperature, envelope mass, and to- 
tal mass. 9 Since the shock brea kout and plateau de- 
pend on M e j, E, and i? pr e SN (e.g.. lEastman et ailll994l : 
IMatzner fc McKee I H999) , our calculation achieves self- 
consistent multicolor LCs from the shock breakout to 
plateau and tail. In this Letter we present a SN explosion 
of a star with a zero-age main-sequence mass Mzams = 
2OM having a presupernova mass M prc sN = 18.4M©, 
H envelope mass M cnv = 13.4M0, and presupernova ra- 
dius -RprcSN = 800-R©. An extensive investigation will be 
presented in a forthcoming paper. 

3. COMPARISONS WITH OBSERVATIONS 

ISchawinski et all (|2008D and lGezari et all (|2008D found 
a UV brightening at the SNLS-04D2dc position which 
lasts several days from ~ 15 days before the first SNLS 
observation. Although the optical observation of the 
shock breakout is not available, the UV-optical LCs of 
the shock breakout and plateau can be compared with 
the synthetic multicolor LCs. We assume the date of 
shock breakout to be 2453062.2 JD (t = 0) and compare 
the model and observations with reference to the date. 
In this Letter the epochs are described in the observer 
frame. 

The multigroup spherical radiation hydrodynamics cal- 
culation provides wavelength- and time-dependent fluxes 
at the SN surface. When a SN is observed from a given 

9 We note that the presupernova progenitor structure depends on 
the treatments of phys ics, e.g., rotation, mass l oss, mixing length, 
and overshooting (e.g., Limongi & ChicfS 2006). 
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Fig. 1. — Top: Sensitivity curves of multicolor bands (red: FUV- 
band, green: NUV-band, blue: g-band, magenta: r-band, cyan: 
i-band, and black: 2-band). For illustration purpose, each filter 
band is blueshifted to the rest frame to compensate for z = 0.185. 
Bottom: Evolution of intrinsic SEDs of a SN IIP model with 
Mzams = 2OM and E = 1.2 X 10 51 ergs at t = day (red), 
0.5 days (green), 2 days (blue), 10 days (magenta), 20 days (cyan), 
50 days (black), and 90 day s (orange). A synt hetic non-LTE spec- 
trum is also shown (violet. fGezari et al.H2008T ). The inset enlarges 
the UV emission at t = 0, 0.5, and 2 days and the non-LTE spec- 
trum. 

direction, lights from different parts of the SN surface are 
radiated at different time and at different radii (for de- 
tails. see Klein fc Chevalierlll978l: llmshennik et al.lll98"il 
lEnsman fc Burrows! 119921 : iBlinnikov et alJ 120021 12003) . 
Thus, we take into account a light travel time cor- 
rection and limb darkening in t he Eddington approxi- 
mation (|Klein fc Chevalierl 119780 . Figure Q] shows cor- 
rected wavelength-dependent luminosities L\ for a SN 
IIP model with M ZA ms = 20M Q and E = 1.2 x 10 51 ergs. 

The SED at t = 2 days i s compared with synthetic 
non-LTE spectrum (55 .6 hr, iGezari et al.l 120081 : see also 
iDessart fc Hilli cr 2005) which gives similar optical color. 
Although the epochs are different because of adopting 
different progenitor models, the UV SED and spectrum 
derived from the independent calculations are distinctly 
consistent. The consistency justifies both theoretical cal- 
culations. 

In order to predict multicolor observations from the 
multicolor theoretical model, the model lights are di- 
luted, 10 redshifted, reddened, and then convolved with 
the sensitivities of th e satel lite and telescope (GALEX: 
iMorrissev et alJ 120051 I2007L the MegaPrime/MegaCam 
on th e Canada-France-H awaii Telescope (CFHT) for 
SNLS: lAstier et al.l [2006h . For illustration purpose, the 
sensitivity curves blueshifted to the rest frame to com- 
pensate for z — 0.185 are shown in the top panel of Fig- 
ure [TJ In this Letter the bands are described in the ob- 
server frame. 

3.1. Ultraviolet light curves of shock breakout 

10 The distance is derived with the five-year resu lt of Wilkinson 
Microwave Anisotropy Probe (Komatsu ct al. 2009). 
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Fig. 2. — Comparison between the NUV (top) and FUV (bot tom) 
observations (points, SNLS-04D2dc, Schawinski ct al. 2008|) and 
the SN IIP model without the host galaxy extinction (dot- dashed 
line) and reddened for the host galaxy extinction with i?s_v host = 
0.06 mag (dashed line), 0.14 mag (solid line), and 0.22 mag (dotted 
line) . 

The UV light allows for direct observations of the 
shock breakout but is strongly reduced by extinction. 
Thus it is crucial to estimate correctly the extinction 
in the host galaxy and our Galaxy. While the color 
excess of our Galaxy Eb-v,Gh\ is taken from Schlegel 
et al. (1998, E B ~v,G&\ = 0.02 mag), it is difficult to es- 
timate the color excess and thu s extinction of the host 
galaxy. [Schawinski et all (|2008[ ) estimates the color ex- 
cess of the host galaxy £?B-v,ho8t & t the SN location 
from the Balmer decrement as EB-v,host = 0.14 mag. 
However, they caution that the uncertainty of the to- 
tal extinction is as much as a factor of two and fur- 
ther note t hat the estimate fro m the empirical relation 
of SNe IIP (jNueent et all 120061) is consistent with both 
of E B -v,host = 0.14 mag and 0. 

Assuming the Small Mage llanic Clo ud (SMC) redden- 
ing law for the host galaxy (Pci 1992) and Es-v,host = 
0.14 mag, the total extinction at the effective wave- 
lengths of the far and near UV (FUV and NUV) filters 
of the GALEX satellite are as large as Apuv = 2.38 mag 
and Anuv = 1-51 mag, respectively. Hereafter, we call 
these values "standard" ex tinction. Althoug h the v alues 
are slightly different from iSchawinski et al.l (|2008l ). the 
total extinction integrated over each band depends on 
the intrinsic spectrum and varies with time as the spec- 
trum changes. The variations of extinction in the UV 
bands are relatively large; for example, ~ 0.3 mag in the 
FUV band and ~ 0.1 mag in the NUV band from t = 
to 20 days for the SN IIP model with A/ Z ams = 20M Q 
and E = 1.2 x 10 51 ergs. 

Because of the large uncertainty, we assume sev- 
eral values for the color excess of the host galaxy as 
follows: Eb-v host = mag referring to the case 
of no extinction in the host galaxy, EB—VJaast = 
0.06 mag giving half of the standard extinction in the 
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Fig. 3. — Bolometric LC (top), color temperature evolution (mid- 
dle), and photospheric velocity evolution (bottom) of the SN IIP 
model (lines). The insets in the top and middle panels enlarge the 
phase of shock breakout. 

NUV band, -Eg-y^ost = 0.14 mag being the stan- 
dard extinction, and Es-v.host — 0.22 mag giving 
double of the standard extinction in the NUV band, 
which lead to (Anuv,^4fuv) = (0.18 mag, 0.15 mag), 
(0.75 mag, 1.10 mag), (1.51 mag, 2.38 mag), and 
(2.27 mag, 3.65 mag), respectively. Here, we assume 
the SMC reddening law for the host galaxy. 

Figure [2] shows comparisons of UV LCs with the model 
with M ZAMS = 2OM and E — 1.2 x 10 51 ergs. The 
model LCs are consistent with the observations within 
the uncertainty, while they are slightly fainter than the 
observations for Eb -v,host = 0.14 mag. The explosion 
energy of SNLS-04D2dc is consistent with the canonical 
value of the explosion energies of core-co llapse SNe [e.g., 
SN 1987A: E = (1.1 ± 0.3) x 10 51 ergs, iBlinnikov et~aTl 
I2000j . Although the 56 Ni- 56 Co radioactive decay does 
not contribute to the shock breakout, we expediently as- 
sume a canonical 56 Ni ejection without mixing to the 
envelope [the ej ected 56 Ni mass M( 56 Ni) = 0.07M©, 
e.g., SN 1987A: IBlinnikov et al.ll2000j . and thus M oj = 
16.9M to yield O.O7M of 56 Ni. 

The second peak in the NUV LC at t ~ 3 days is repro- 
duced by the model and explained by t he shift of the peak 
wavele ngth as [Schawins ki et ail (|2008f ) and lGezari et al.l 
( 2008) suggested. The bolometric LC and the evolution 
of color temperature are shown in Figure [3] Figure [3] 
also shows the velocity evolution of photosphere defined 
as a position where the radiation and gas are decoupled. 
Although the bolometric luminosity declines monotoni- 
cally after the shock breakout, the radiation energy in 
the NUV band increases with time because the peak 
wavelength shifts long (Fig. [J). After the NUV second 
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Fig. 4. — C omparison between the SNLS observations (points, 
SNLS-04D2dc, ISchawinski et alj |2008l and the SN IIP model red- 
dened for the host galaxy extinction with Es-v,host = 0.14 mag 
(lines) (red: (/-band, green: r-band, blue: i-band, magenta: 2-band 
in AB magnitude system). 

peak, the color temperature decreases further and the 
peak wavelength shifts to the optical bands. The shift is 
caused by not only the decreasing temperature of the SN 
ejecta but also an enhancement of the metal absorption 
lines due to the low temperature. As a result, the UV 
luminosity declines monotonically after the NUV second 
peak. The model also predicts a second peak in the FUV 
band but the brightening is obscured because of the low 
signal-to-noise ratio. 

3.2. Optical light curves at plateau stage 

Thanks to the multigroup radiation hydrodynamics 
calculations, subsequent evolutions of multicolor lights 
are obtained and compared with the SNLS optical ob- 
servations. Figure |4] shows comparisons of the g-, r-, i-, 
and z-band LCs. Here, we adopt EB-v,has% = 0.14 mag 
and the SMC reddening law for the host galaxy. The 
total extinction at the effective wavelengths of the g-, r- 
, i-, and z-band filters of the MegaPrime/MegaCam on 
CFHT are as large as A g = 0.64 mag, A r = 0.47 mag, 
Ai = 0.36 mag, and A z = 0.28 mag, respectively. 

As SED peaks in the g-band at t ~ 20 days (Fig. [T]) , 
the g-band LC peaks at t ~ 20 days. After this epoch, 
the SED becomes red with time and the blue edge of the 
SED enters in the g-band. Thus, the g-band LC declines 
more rapidly than the other optical-band LCs. On the 
other hand, the r-band LC declines gradually and the 
decline rate slightly changes at t ~ 55 days when the 
SED peak enters in the r-band. The i- and z-band LCs 
brighten by t ~ 60 days due to the shift of the peak 
wavelength. 

Although the overall multicolor LCs are well repro- 
duced by the model, the z-band model LC is brighter 
than the observation at t ~ 90 days. The discrepancy 
would be improved if we add more physics to STELLA. 
STE LLA currently has ar ound 1.5 x 10 5 spectra l lines 
from lKurucz fe Belli (|1995l ) and lVerner et all (|19960 and 
a rather poor line list in near infrared, while tens of 
millions of spectral lines from large Kurucz lists are be- 
ing included (E. Sorokina, private communication). The 
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Fig. 5. — Apparent B-band light curve of a shock breakout in 
AB magnitude system for a SN being identical to SNLS-04D2dc 
at z = 1. We assume Bs_v,Gal = 0-02 mag, -Es-^host = 
(solid line) and 0.14 mag (dashed line), and the SMC red- 
dening law for the host galaxy. A limiting magnitude for a 
4<t de tection in a hour of integration with SUBARU/Suprime- 
Cam (Mivazaki ct al. 2002) is also shown (dotted line), which 
is calculated with Subaru Imaging Exposure Time Calculator 
(http://www.naoj.org/cgi-bin/spcam_tmp.cgi) assuming 0.5" see- 
ing, 1.5" aperture, and '6 days trom New Moon. 

expansion of line lists, together with taking into ac- 
count non-equilibrium effects in an important coolant 
like Ca II, m ay influence the z-band flux appreciably 
(|Kasenll2006D . 

No mixing of 56 Ni leads to no contribution from the 
radioactive heating to the LC at t < 100 days, while 
the radioactive decay could power the LC from the early 
phases if 56 Ni is mixed to the envelope. Because of the 
lack of observations at t > 100 days, M( 56 Ni) of SNLS- 
04D2dc cannot be constrained. We note that the avail- 
able observations can be reproduced equally well even if 
no 56 Ni is ejected. 

4. CONCLUSIONS & DISCUSSION 

We present a multicolor LC model with the multi- 
group radiation hydrodynamical code STELLA. The 
model reproduces the multicolor UV-optical LCs of the 
shock breakout and plateau of SNLS-04D2dc consis- 
tently. Since the shock breakout is more sensitive to 
E and -RnrsSN than the plateau ((Eastman et al.1 Il994t 
iMatzner fc McKee1ll999f ). the SN and progenitor prop- 
erties are constrained more tightly than only with the 
plateau observation. We clarify the properties of SNLS- 
04D2dc: the progenitor is a star with Mzams = 20M Q , 
MprcSN = 18.4M©, and i? proS N = 800i? Q , and the SN 
has the canonical explosion energy E — 1.2 x 10 51 ergs. 
Because of no observations at t > 100 days, we cannot 
constrain M( 56 Ni) of SNLS-04D2dc, while we expedi- 
ently adopt M e j = 16.9M Q and M( 56 Ni) = 0.07M Q . 

The second peak of the NUV LC is reproduced by the 
shift of the peak wavelength, in spite of the monotonic de- 
cline of the bolometric LC. Although the signal-to-noise 
ratio of the FUV LC is quite low, the model is also in 
an agreement with the FUV LC. The consistency sup- 
ports the origin of the NUV second peak. When only the 
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monotonic light curve is obtained, only two observational 
quantities, i.e., flux and duration, are available and the 
three properties characterizing the shock breakout, M e j, 
-RpreSN, and E, are degenerated. The multicolor obser- 
vation enables us to constrain the color temperature and 
thus to resolve the degeneracy. In order to clarify the SN 
properties in detail, it is required to employ the multi- 
group radiation hydrodynamical calculations because the 
spectra of the shock breakout slightly deviate from the 
blackbody. 

The successful reproduction of UV-optical LCs of 
SNLS-04D2dc and the consistency with the non-LTE 
spectral calculation justify the multigroup radiation hy- 
drodynamics calculation with STELLA and indicate 
that the calculations are capable of predicting the multi- 
color observations of the shock breakout and plateau of 
SNe IIP. If a SN being identical to SNLS-04D2dc takes 
place at z = 1 in the direction with Eb-v,g&\ = 0-02 mag, 
the SN brightness will be me ~ 26.4 mag without the 
host galaxy extinction or me ~ 27.6 mag with the host 
galaxy extinction with Es-v,host = 0.14 mag and the 
SMC reddening law (Fig. [5]). It is bright enough to be 
detected with 8m-class telescopes. The result envisages 
that future deep and wide surveys by, e.g., SkyMapper, 
Panoramic Survey Telescope and Rapid Response Sys- 
tem (Pan-STARRS), Large Synoptic Survey Telescope 
(LSST), and SUBARU/Hyper Suprime-Camera (HSC), 
will find a large number of shock breakouts. The large 
sample will import totally-new knowledge about cosmic 
evolution histories, e.g., a core-collapse SN rate and a 
star formation rate, at z > 1. Therefore, although the 
detections of Type Iln SNe at z > 2 are recently reported 
(jCooke et al. 2009), we still emphasize a great potential 



of shock breakouts for direct detection of the most com- 
mon core-collapse SNe, SNe IIP, at z > 1. 

Finally, we point out that the estimate of the host 
galaxy extinction is crucial to estimate the intrinsic lu- 
minosity of shock breakout. This is because the shock 
breakout of SN IIP emits the radiation energy mainly in 
UV easily reduced by the interstellar extinction. With 
the use of the sensitiveness, if the observations of shock 
breakout, plateau, and tail are available with high signal- 
to-noise ratios, finding a set of a SN model and ex- 
tinction, which consistently reproduces the observations, 
might provide a new constraint on the host galaxy ex- 
tinction. 
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